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ABS'^ACT 

~^The  physicochemical  phenomena  associated  with 
fluid  flow  In  an  evaporating  ultrathin  film  of  an 
Ideal  binary  mixture  In  the  contact  line  region  are 
discussed.  The  effect  of  composition  and  temperature 
gradients  on  surface  shear  1$  analyzed  and  the 
results  are  compared  with  experimental  trends. 
Although  simple  In  concept,  the  constant  vapor 
pressure  boundary  condition  appears  to  be  very  useful 
In  delineating  the  relative  effects  of  surfar^  ten¬ 
sion,  temperature  and  composition  on  fluid  flow  due 
to  surface  tension  gradients  In  general, 

NOHENCLATURE 


'slv 


■  Hamaker  constant 

■  -*slv'®’ 

•  parameter  In  Equation  (3) 

•  Intermolecular  distance 


g  ■  acceleration  due  to  gravity 

K  •  curvature 

n  •  parameter  In  Equation  (3) 


V  •  velocity  in  the  y  direction 

X  ■  distance  parallel  to  the  main  direction  of 

flow 

X  >  concentration  of  component  •  (n  ■  1  or 

"  n  -  2) 

y  ■  distance  perpendicular  to  the  main  direc¬ 

tion  of  flow 

Greek 

r  ■  mass  flow  rate 

<  'film  thickness 

e  ■  angle  of  Inclination 

II  •  absolute  viscosity 

V  ■  kinematic  viscosity 

p  ■  density 

a  >  surface  tension 

T  ■  Shear  stress  In  x  direction  on  a  surface 

perpendicular  to  the  y  direction 


p  •  pressure 

p^  •  pressure  In  the  vapor  space 

p^  •  vapor  pressure 

p*  >  saturation  vapor  pressure 

P  •  p  ♦  a 

T  >  absolute  temperature 

u  ■  velocity  In  the  x  direction 


f  ■  potential  energy  function  per  unit  volume 

•g  ■  potential  energy  per  unit  volume  of  a  sur¬ 

face  molecule  of  a  bulk  phase 

Subscripts 

1  ■  lighter  component  In  a  binary  mixture 

2  ■  heavier  component  In  a  binary  mixture 

Is  ■  liquid-solid 


n 


liquid 

solid 


ss 

siv  ■  solid-liquid-vapor 

INTRODUCTION 

In  a  recent  set  of  publications,  the  use  of  a 
scanning  microphotometer  to  help  determine  the  heat 
transfer  characteristics  of  a  thin  stationary  eva¬ 
porating  film  of  liquid  in  the  film  thickness  range  i 
<10~^  m  was  described  [1-5].  A  simple  diagram  of  the 
system  is  presented  in  Fig.  (1).  The  microphotometer 
was  used  to  determine  the  thickness  profile  by 
measuring  the  relative  location  of  Interference 
fringes  formed  in  the  film  by  the  reflection  of 
light.  Since  various  mechanisms  causing  fluid  flow 
are  a  function  of  the  film  profile  it  can  be  used  to 
evaluate  their  effect  on  evaporative  heat  transfer. 
Mixtures  of  alkanes  of  various  weight  percent  were 
studied.  The  initial  studies  were  concerned  with 
relatively  pure  systems  and  were  designed  to  measure 
the  effect  of  evaporation  rate  on  the  thin  film  pro¬ 
file  [1].  Fluid  flow  resulting  from  a  capillary 
pressure  gradient  (due  to  a  curvature  gradient)  and 
from  a  disjoining  pressure  gradient  (due  to  a 
thickness  gradient)  were  identified.  The  evaporating 
thin  films  were  found  to  be  very  stable  with  fluid 
flowing  towards  the  heat  source.  In  addition,  small 
percentages  of  a  second  component  were  experimentally 
found  to  be  very  important.  It  became  apparent  that 
even  small  quantities  of  an  Impurity  with  a  substan¬ 
tially  lower  vapor  pressure  would  concentrate  as  a 
result  of  distillation  as  the  fluid  flowed  towards 
the  heat  source.  Subsequent  studies  demonstrated 
that  small  changes  in  the  bulk  composition  signifi¬ 
cantly  altered  the  characteristics  of  the  transport 
processes  in  the  contact  line  region  of  evaporating 
thin  films  [2-5].  The  curvature  gradient  at  the 
liquid-vapor  interface  was  found  to  be  a  strong  func¬ 
tion  of  evaporation  rate  and  bulk  composition,  and  to 
be  strongly  coupled  to  the  concentration  gradient. 
Concentration  and  temperature  gradients  apparently 
caused  interfacial  shear  stresses  and  fluid  flow  pat¬ 
terns  that  enhanced  contact  line  stability.  Although 
the  thickness  profile  could  be  easily  measured  opti¬ 
cally.  the  concentration  and  temperature  profiles  in 
a  region  less  than  10'^  m  in  length  could  not  be 
easily  obtained  experimentally.  Therefore,  extensive 
modeling  of  the  fluid  flow  mechanisms  dependent  on 
these  variables  is  needed.  Herein,  a  simple  model 
that  allows  the  relative  importance  of  various  physi¬ 
cochemical  phenomena  on  fluid  flow  in  an  evaporating 
thin  film  to  be  evaluated  is  presented.  The  effect 
of  composition  and  temperature  gradients  on  the  sur¬ 
face  shear  is  emphasized.  The  experimental  results 
on  evaporating  thin  films  of  mixtures  [2-5]  parallels 
in  many  ways  the  pioneering  experimental  results 
obtained  by  Bascom,  Cottington  and  Singleterry  on 
spontaneous  spreading  [6].  This  research  on 
spreading  concerned  extremely  low  evaporation  rates. 
However,  surface  tension  gradients  due  to  con¬ 
centration  gradients  were  Identified  as  the  driving 
force.  This  work  was  extended  by  Ludviksson  and 
Lightfoot  who  used  a  nonisothermal  plate  with  its 
warm  end  partially  immersed  and  fluid  flowing  away 
from  the  heat  source  [7].  Me  note  that  there  are 
many  technologies  in  heat  transfer  that  rely  on  eva¬ 
porating  thin  films  that  are  related  to  the  basic 
problem  addressed  here. 


Fig.  1.  Evaporating  Extended  Meniscus  on  Inclined 
Flat  Plate  Immersed  in  a  Pool  of  Liquid. 

FLUID  FLOW  IN  AN  EVAPORATING  FILM 

A  continuum  model  is  used  to  discuss  the  physi¬ 
cochemical  phenomena  of  importance  to  fluid  flow  in  a 
thin  film  with  a  varying  thickness  even  though  the 
thickness  in  a  portion  of  the  flow  field  can  be  less 
than  the  minimum  thickness  where  the  continuum  mdel 
is  obviously  justified.  However,  we  feel  that,  as  a 
minimum,  the  resulting  insights  Justify  this  presump¬ 
tive  extension  to  include  ultra-thin  regions  in  the 
initial  discussion.  Me  also  note  that  the  second 
order  effect  of  diffusion  is  neglected  in  this  paper 
even  though  the  effect  of  concentration  on  the  sur¬ 
face  tension  gradient  is  discussed.  Neglecting  iner¬ 
tia  terms  and  the  vertical  component  of  velocity  the 
following  simplified  Navier-Stokes  equation  can  be 
written  for  the  velocity  in  the  evaporating  thin  film 
shown  in  Fig.  (1). 

In  which 

♦  •  -pgx  sin  e  ♦  eg  -  .Sjj-  (2) 
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Herein  we  use  a  potential  energy  function  per  unit 
volume  which  is  a  function  of  the  film  thickness,  y  ■ 
<(x).  to  represent  the  difference  in  behavior  of  a 
thin  film  relative  to  that  of  a  semi-infinite  bulk 
phase  [8-13].  eg  t***  potential  energy  per  unit 
volume  of  a  surface  molecule  of  the  bulk  phase. 
Although  the  results  are  the  same  for  all  practical 
purposes  as  those  obtained  using  a  pressure  Jump  at 
the  interface,  e.g.  [14,15],  this  approach  circum¬ 
vents  the  need  to  discuss  disjoining  pressure  expli- 


COMPARISON  OF  EFFECTS  DUE  TO  SURFACE  AND  BULK 
POTENTIAL  GRADIENTS 


citly.  The  Hamaker  constant  for  a  spreading  system 
Is  Aj^y  <  0.  Equation  (1)  can  be  solved  for  the 
velocity  distribution  using  the  following  boundary 
conditions. 

y  •  0,  u  «  0  (5) 


y  •  «(x).  ^yx  ■  0  • 


y  •«(x).P.p„  +  *8-^-0  K  (7) 

& 

The  letter  x  without  a  subscript  refers  to  position, 
whereas,  the  letter  x  with  a  subscript  refers  to  con¬ 
centration.  The  second  boundary  condition  can  have  a 
large  effect  on  the  velocity  distribution  because 
large  gradients  in  the  concentration  and  temperature 
are  possible  In  evaporating  thin  films.  Even  a  bulk 
fluid  that  Is  relatively  pure  can  develop  large  con¬ 
centration  gradients  In  an  evaporating  thin  film  due 
to  distillation.  For  extremely  thin  films  the  sur¬ 
face  tension  Is  also  a  function  of  the  film 
thickness.  These  additional  effects  due  to  the 
effect  of  the  film  thickness  on  the  surface  tension 
are  beyond  the  scope  of  the  present  paper.  The  total 
pressure  In  the  vapor  space,  Po>  equal  to  the 

vapor  pressure  because  It  Includes  a  variable  non¬ 
condensible  component.  We  note  that  the  dispersion 
term  gives  a  bulk  effect  and  the  surface  tension  a 
surface  effect  on  fluid  flow  In  the  thin  film.  Using 
Equations  (1-7)  we  obtain  Equation  (8)  for  the 
Integration  of  pu  across  the  film. 

r  ■  ^  [1.54 *^o'  +  +  K<j'  ♦  oK'  -  nB4‘*"*^'«' 

(I)  (H)  (III)  (IV)  (V) 

-  pg6'  cos  e-^pg  sin  e]  (8) 


The  following  model  of  Equation  (9)  Is  for¬ 
mulated  to  evaluate  the  significance  and  direction  of 
surface  tension  driven  flow  for  a  qualitative 
understanding  of  past  experimental  data.  Since  the 
surface  tension  Is  a  function  of  both  the  con¬ 
centration  and  the  temperature  It  Is  advantageous  to 
use  an  equilibrium  boundary  condition  of  a  constant 
vapor  pressure  for  the  most  volatile  component. 
Physically  It  seems  reasonable  to  assume  that  the 
vapor  pressure  for  a  short  distance  along  the  vapor- 
liquid  Interface  would  be  constant  and  approximately 
equal  to  the  equilibrium  vapor  pressure  If  the  Inter- 
facial  resistance  to  evaporation  Is  relatively  small 
and  If  tangential  vapor  flow  occurs  rapidly.  In 
Appendix  I,  a  boundary  condition  based  on  effusion  Is 
found  to  be  approximately  equal  to  a  constant  vapor 
pressure  boundary  condition  for  the  conditions  under 
study.  We  note  that  the  constant  vapor  pressure 
boundary  condition  Is  not  a  constant  evaporative  heat 
flux  boundary  condition  because  the  variation  of  com¬ 
position  across  the  film  Is  not  addressed. 

Therefore,  the  concept  appears  to  be  of  general  use 
when  discussing  the  surface  shear  direction.  In 
addition,  the  assumption  that  the  second  component  of 
the  binary  mixture  being  modeled  has  a  negligible 
vapor  pressure  Is  In  agreement  with  some  past  experi¬ 
mental  systems  [2-5], 

The  variation  with  position  of  the  surface  ten¬ 
sion  of  a  binary  mixture  can  be  written  as 
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In  which  the  letter  x  without  subscript  refers  to 
position. 

Using  Equation  (13)  In  this  equation  gives 
Equation  (14)  for  a  binary  mixture. 


(VI)  (VII) 

The  prime  refers  to  differentiation  with  respect  to 
X.  In  this  paper  we  will  focus  on  the  first  three 
terms  In  the  above  brackets  which  account  for  flow 
due  to  gradients  In  the  Intermolecular  forces 
resulting  from  gradients  In  the  concentration  and 


temperature: 

(I)  •  1.5  4"'  o' 

(9) 

(II)  -  4'"  B' 

(10) 

(III)  ■  Ko* 

(11) 

The  fourth  and  fifth  terms  primarily  depend  on  the 
film  profile  which  can  be  experimentally  determined 
and  have  been  found  to  be  functions  of  the  com¬ 
position  and  evaporation  rate  [1-5].  The  final  two 
terms  are  relatively  small  except  when  r  *  0  and  can 
be  evaluated.  Experimental  data  have  demonstrated 
that  the  third  term  Is  small  compared  to  the  first 
term  since  K  Is  usually  small  compared  to  1.5  4*1 
[1-5].  At  present,  the  effects  of  the  first  two 
terms  on  fluid  flow  are  the  primary  unknowns  and  will 
be  evaluated  below. 


o  ■  XjOj(T)  ♦ 

3x  ^*1  ?r  *  *2  3x  *  <"i  •  "2^  W 

Using  Equation  (15),  which  Is  an  extension  of 
Raoult's  Law  that  Includes  the  effect  of  a  "pressure 
Jump"  at  the  HqUld-vapor  Interface,  aP  ■  Po*^» 
the  variation  of  vapor  pressure  with  temperature, 
pressure  and  composition  gives  Equation  (16)  for  the 
variation  of  vapor  pressure  with  position.  The 
current  model  applies  to  that  portion  of  the  thin 
film  where  the  curvature  Is  suficlently  small  and  the 
thickness  Is  sufficiently  thick  so  that  the  effects 
of  curvature  and  the  Interfaclal  force  due  to  the 
solid  on  the  vapor  pressure  can  be  neglected. 
Therefore,  the  "pressure  Jump"  at  the  Interface  Is 
taken  to  be  zero. 

V 

Pvi  ‘  *i'*5i  **'* 


for  aP  *  0 


(16) 


Imposing  a  constant  component  vapor  pressure  boundary 
condition  gives  Equation  (17)  which  can  be  used  along 
with  Equation  (14)  to  give  Equation  (18)  for  the 
variation  of  surface  tension  with  position. 


dx 
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The  evaluation  of  the  various  terms  In  Equation  (18) 
Is  discussed  In  Appendix  II  and  the  results  are  pre¬ 
sented  In  Figures  (2-4).  The  particular  systems 
selected  were  studied  In  Refs.  [2-5].  Before 
discussing  these  results,  the  relative  Importance  of 
the  bulk  flow  term,  Eq.  (10),  will  be  analyzed  In  the 
next  paragraph. 


Using  Equation  (19)  for  the  substrate  portion  of 
the  liquid-solid  Hamaker  constant  [16]  In  Equation 
(3)  gives  Equation  (20). 

A„  .  (W 

'.I.  ■  *n  -  <“l 

Taking  the  derivative  of  Equation  (20)  with  respect 
to  position  while  assuming  that  Ajj  Is  not  a  function 
of  position  gives  Equation  (21). 


Tx —  ‘  -3x-  J 


(21) 
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Using  A..  ■  24ii0  «  for  the  liquid  film  portion  of  the 
Hamaker  Constant  [17]  In  Equation  (21)  gives 


.  24,0^  [1  -  ||  (22) 

The  distance  0  Is  the  order  of  intermolecular  distan¬ 
ces.  Taking  the  ratio  of  Equations  (9)  and  (10) 
using  Equations  (3)  and  (22)  gives  Equation  (23) 
which  demonstrates  that  the  effect  due  to  a  variation 
of  surface  tension.  Equation  (9),  Is  usually  more 
significant  than  the  effect  due  to  a  variation  In  the 
Hamaker  constant.  Equation  (10),  since  S  »  0. 

He  note  that  the  additional  effect  of  thickness  on 
the  vapor  pressure  In  Equation  (16)  still  needs  to  be 
evaluated  for  completeness. 

OISUCSSION  OF  RESULTS 


flow  changes  sign  In  a  distilling  thin  film  as  It 
flows  towards  the  heat  source  along  a  constant  vapor 
pressure  line  with  (02  >  o,)*  Assuming  that  the  flow 
starts  In  a  region  where  surface  shear  enhances  the 
flow  towards  the  heat  source,  do/dT  >  0,  a  com¬ 
position  Is  reached  In  the  flowing  film  at  which  flow 
reversal  can  occur  In  the  film.  Assuming  further 
that  there  Is  an  Initial  transient  period  during 
which  the  flow  at  this  point  Is  greater  than  the  eva¬ 
poration  rate  above  It,  a  liquid  reservoir  In  the 
shape  of  a  bulge  should  form  since  some  evaporation 
In  this  region  would  lead  to  a  concentration  at  a 
higher  temperature  closer  to  the  heat  source  at  which 
the  surface  shear  opposes  the  flow.  This  result 
appears  to  be  of  general  Importance  to  the  formation 
and  movement  of  lenses.  The  results  presented  In 
[2-5]  and  In  the  figures  demonstrate  that  the  for¬ 
mation  of  a  bulge  Is  a  function  of  the  bulk  con¬ 
centration,  temperature  level  and  ao:  an  Increase  In 
the  concentration  of  a  lower  boiler  with  a  signifi¬ 
cantly  different  surface  tension  In  the  bulk  leads  to 
an  Increase  In  the  propensity  to  form  a  pool  above 
the  meniscus  for  a  limited  range  of  the  surface  tem¬ 
perature  and  concentration.  He  note  that  the 
coupling  between  the  initial  composition,  the  frac¬ 
tion  distilled  and  the  final  composition  Is  not  pre¬ 
sented  In  the  figures. 


The  model  predicts  the  following  experimentally  Fig.  2.  Surface  Tension  Derivative  (Coefficient  In 

observed  trends  [2-5]:  Equation  (18))  Versus  Temperature  for 

Hexane-Octane  Nixtures. 

The  Formation  of  a  Bulge  In  the  Thin  Film  Portion  of 
I  the  Extended  Meniscus. 


The  results  presented  In  Figures  (2)  and  (3) 
demonstrate  that  the  effect  of  surface  shear  on  fluid 


X,  — 1.0 


3:  Pvi  *  292  mmHg 
2:  Pvt  “  72  mmHg 
1:  Pvt  *  11  mmHg 


Fig.  3.  Surface  Tension  Derivative  (Coefficient  In 
Equation  U8))  Versus  Temperature  for 
Oecane-Tetradecane  Mixtures. 


The  Formation  of  a  Bulge  would  not  be  Observed  with  a 
Very  Pure  Fluid. 

Since  distillation  Mould  change  the  composition 
of  a  very  pure  fluid  very  slowly,  this  case  Is  repre¬ 
sented  by  XX  ■  1.0  In  the  figures.  The  results  pre¬ 
sented  In  Figure  (4)  demonstrate  that  the  surface 
shear  Is  away  from  the  higher  temperature  If  the  dif¬ 
ference  In  the  surface  tensions,  ia  •  (o,  -  o.)* 
relatively  small,  which  is  the  usual  case  for  relati¬ 
vely  pure  samples.  In  these  systems  flow  towards  the 
heat  source  Is  due  to  Terms  (IV  $  V)  In  Equation  (8). 
We  note  that  small  residual  lenses  were  observed  In 
[1]. 

The  Curvature  Gradient  Is  a  Strong  Function  of  the 
Bulk  Composition. 

Equation  (8)  demonstrates  that  the  curvature 
gradient  and  the  surface  tension  gradient  are 
strongly  coupled.  Therefore,  a  change  in  the  surface 
tension  gradient  for  the  same  mass  flow  rate  could 
appear  as  a  change  In  the  curvature  gradient  In  the 
reported  experimental  results.  When  the  experimental 
results  were  analyzed  In  Ref.  [2-5]  using  only  Terms 
(IV  -  VII)  In  Equation  (8),  relatively  more  flow 
apparently  occurred  with  the  pure  system  than  with 
the  mixture  for  the  same  heat  flux  based  on  the  tem¬ 
perature  measurements.  This  apparent  difference 
could  be  the  result  of  the  coupling  between  the 
various  terms  In  Eq.  (8). 


Fig.  4.  x^  &o  Versus  Temperature  for  o'  ■  0  and 

do.  do*  . 

TT  "  -  .lO*'  H/mK  for  Hexane. 

CONCLUSIONS 

A  simple  model  based  on  continuum  mechanics, 
interfaclal  phenomena,  a  constant  vapor  pressure 
boundary  condition,  and  an  extension  of  Raoult's  Law 
qualitatively  describes  the  following  trends  experi¬ 
mentally  observed  In  experiments  of  fluid  flow  In 
evaporating  thin  films  In  the  contact  line  region: 

1)  The  effect  of  surface  shear  can  change  sign  In  a 
distilling  thin  film  as  It  flows  towards  the 
heat  source  along  a  constant  vapor  pressure 
line; 

2)  There  Is  a  significant  range  of  composition  and 
temperature  over  which  surface  shear  enhances 
flow  towards  the  heat  source; 

3)  Surface  shear  Inihibits  flow  towards  the  heat 
source  In  a  relatively  pure  fluid  when  the  dif¬ 
ference  In  component  surface  tensions  Is  small; 
and 

4)  The  surface  tension  gradient  and  the  curvature 
gradient  are  strongly  coupled. 
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APPENDIX  I 

Schrage  [18]  used  the  concept  of  effusion  to 
derive  an  equation  for  the  net  mass  flux  of  matter 
crossing  a  liquid-vapor  Interface  which  has  been  used 
extensively,  e.g.,  [13,19-21].  The  mass  flux  due  to 
effusion  at  an  Interface  Is  proportional  to  py/ZT. 
Taking  this  to  be  constant  at  the  liquid-vapor  Inter¬ 
face  an  alternate  equation  to  Equation  (18)  can  be 
derived  for  the  variation  of  surface  tension  with 
position. 
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In  place  of  Eq.  (15)  gives 


1  .  X  r  1  1  1  fA  2 

w  *1  ■^TT-  •  zr-J 

For  temperatures  of  Interest  herein,  the  second  term 
in  the  brackets  In  Eq.  (A-2)  Is  small  compared  to 
the  first  term, 

APPENDIX  II 

In  order  to  evaluate  the  surface  tension  as  a 
function  of  temperature  the  following  equations 
recommended  by  Jasper  [22]  were  used  (Temperatures, 
T^,  are  In  Celsius  units): 

Hexane:  o  •  (20.44  -  0.1022  T^)  x  10'^  N/m 

Octane:  a  -  (23.52  -  0.09509  T^)  x  10‘^  N/m 

Oecane:  o  •  (25.67  -  0.09197  T^)  x  10'^  N/m 

Tetradecane:  a  «  (28.30  -  0.08682  T^)  x  10“^  N/m 

The  following  Antoine  equation  was  used  to 
represent  the  saturation  vapor  pressure  In  imiHg  for 
the  alkanes: 

In  p*  •  A  -  (A-3 

Differentiation  of  Eq.  (A-3}  leads  to 


The  following  Antoine  constants  were  found  In  the 
literature: 


Decane  [231  I  16.0114  |  3456.80 


Hexane  [23 
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